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ABSTRACT 


Context. Albireo is a well-known bright visual double star. It is still unclear if the components A and B form a gravitationally bound 
system. The component Albireo A is itself a binary star. The orbital parameters of the Albireo Aa, Ac system have been determined 
only recently. Thus, Albireo is still of interest for current research. 

Aims. We aim to present evidence for the detection of a new member in the Albireo system. Furthermore, we aim to determine the 
orbital parameters and to find further conclusions for the Albireo system. 

Methods. We used spectroscopic observations of Albireo A obtained with the TIGRE telescope and determined the radial velocities 
during a period of over three years. We analyzed the radial velocity curve with RadVel to determine the orbital parameters. In addition, 
we determined the stellar parameters of Albireo Aa with iSpec. 

Results. We found clear evidence for yet another star in the Albireo system orbiting Albireo Aa with a period of about P = 371 days. 
Several alternative explanations for the periodic radial velocity signal could be discarded. The new companion Albireo Ad is a low 
mass star of about 0.085 Mo. 

Conclusions. We conclude that Albireo is a hierarchical multiple star system and remains an interesting object for future observations 


©ESO 2022 


and studies. 


1. Introduction 


v2 [astro-ph.SR] 14 Mar 2022 


E] Albireo (8 Cyg) is a very bright double star easily visi- 
ble for small-sized telescopes, well-known among amateur as- 
tronomers, and an ideal object for public observations because of 
the very different colors of its components. It has been observed 

* for over more than two centuries (Hass et al.2016). However, 

CO there still exist many open questions about the origin, the compo- 

© nents, and the orbital parameters of the Albireo system. Because 
of this, Albireo is still studied today and has been observed with 


N many different methods (Scardia et al.|2007; 1999 
, z [son et al.|/2013} [Roberts & Mason|[2018 2019). There 


>< is the long-lasting question as to whether the two bright com- 
į ponents Albireo A and B form a gravitationally bound system 
6S or not (Griffin|/1999). The Gaia collaboration 
2016) published improved parallaxes in the recent 
Data Release EDR3 
jet al.|[2021), and it seems that a gravitationally bound orbit is 


still possible because the parallaxes of the components A and B 
lie within a 2o error. It is worth mentioning that Gaia, in gen- 
eral, has difficulties determining parallaxes of bright stars such 
as Albireo (Drimmel et al.!2019). In a thorough spectral analy- 
sis, (2021) have shown that Albireo A and B have 
about the same age and, therefore, might have at least a common 


origin. This conclusion is also supported by the discovery of a 
moving group containing four other fainter stars. 


Key words. Stars: individual: Albireo — binaries: spectroscopic — Techniques: radial velocities 


The star Albireo A is a binary system, which has been known 


for a long time (Maury & Pickering||1897{ 1899). The 


spectrum is a composite of a K3II giant with a B9V compan- 


ion (Markowitz|1969} 1998). This binary system 


has also recently been studied so as to understand its nature and 
to determine its properties (Jack et al.| Bastian & Anton 
[Drimmel et al.2021). (2021) determined 
the orbital parameters of the Albireo Aa, Ac system combining 
a large set of long-term observations of the radial velocity (RV), 
astrometry, and speckle interferometry. They determined the or- 
bital period of the Albireo Aa, Ac system to be about 122 years. 
Because of this long period, the uncertainties of the orbital pa- 
rameters are still quite large. They detected the problem that the 
total mass of the Albireo A system was too high when compared 
to the masses determined by spectroscopic analysis. This "miss- 
ing mass" must be hidden in the star Albireo Ac. They proposed 
that Albireo Ac might be a binary system either consisting of 
two very similar main sequence stars or that there is an invisible 
black hole component. This indicates that there exists at least 
one further object in the Albireo system. 


In addition, there have been reports about the detection of 
another close companion of Albireo Aa using speckle interfer- 


ometry (Bonneau & Foy) 1980 2002). However, the 


detections were only marginal so that the existence of the star 
Albireo Ab has not been confirmed and is still doubtful. 
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Fig. 1. Lomb-Scargle periodogram of the RV measurements of Albireo 
Aa showing a clear peak around 370 days. The graph below shows the 
window function, where the peak disappeared. 


In conclusion, the number of objects in the Albireo sys- 
tem and their properties are still unclear. In this publication, we 
present the discovery of a new low mass star in the Albireo sys- 
tem orbiting Albireo Aa, which we detected with RV measure- 
ments. 


2. Observations of Albireo A 


We obtained a time series of optical spectra with the Heidelberg 
Extended Optical Range Spectrograph (HEROS) echelle spec- 
trograph mounted on the 1.2 m Telescopio Internacional de Gua- 
najuato Robótico Espectroscópico (TIGRE) telescope 
in central Mexico. The spectrograph has a resolution 
of R ~ 20,000 and covers the optical wavelength range from 
about 3800 to 8800 A divided into two channels (red and blue) 
with a small gap of 120 A between the two channels at around 
5800 A. We observed Albireo A for over three years. All spectra 
were obtained with a high signal-to-noise ratio (S/N) between 
200 and 400. One goal of this study was to determine the RVs 
of Albireo Aa. We used the method described in detail in 
(2018). This method has also been successfully used in 
studies of spectroscopic binary stars (Jack et al.[2020} [2022). 

The observed spectra of Albireo A are composite spectra 
that contain contributions from the two stars Albireo Aa and 
Ac (see {Drimmel et al.|(2021) for a complete TIGRE spectrum 
of Albireo A). The contribution from Albireo Ac is only visi- 
ble in the blue channel of HEROS because it is a B type main 
sequence star. The main contribution in the red channel of the 
HEROS spectrograph comes from the K red giant Albireo Aa. 
The method to determine RVs uses only spectral lines in the red 
channel so that all the lines correspond to Albireo Aa. 

We observed Albireo A for over three years and detected a 
clear variation in its RV. We present all values of the RV mea- 
surements in Table|1} where we list the RV in km s~! in terms 
of the Julian date (JD). First, we checked for a possible peri- 
odicity in the RV measurements generating a Lomb-Scargle pe- 
riodogram that we demonstrate in Fig.|1} There is a clear and 
strong peak at around 370 days, which indicates a periodic sig- 
nal in the RV curve. The Fourier analysis gave a best period of 
373.02 days. Below, we plotted the window function, where the 
peak disappeared indicating that it is not just a seasonal effect. 
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Table 1. RV measurements of Albireo Aa in terms of the Julian date 
(JD) obtained with optical spectra using the TIGRE telescope. 


JD RV in km s~! 
2458391.55413 -24.83 +0.13 
2458446.55949 -25.07 + 0.12 
2458527.01935 -25.48 + 0.11 
2458528.01464 -25.56 + 0.11 
2458529.01321 -25.52 + 0.10 
2458577.96934 -25.58 + 0.10 
2458612.92724 —25.47 +0.10 
2458624.88997 -25.38 + 0.11 
2458717.65814 —24.83 + 0.13 
2458748.56362 —24.88 + 0.13 
2458770.59632 —24.95 + 0.13 
2458807.53644 —24.89 + 0.12 
2458898.01421 -25.34 + 0.10 
2458933.92625 —25.44 + 0.10 
2458962.96438 -25.49 + 0.10 
2458991.86885 -25.50 + 0.10 
2459021.80815 -25.38 + 0.10 
2459053.76334 -25.06 + 0.11 
2459054.79683 -24.97 + 0.10 
2459084.75071 —24.75 + 0.14 
2459116.59518 —24.90 + 0.12 
2459145.57801 —24.91 + 0.12 
2459155.55531 —24.89 + 0.13 
2459165.53999 —24.82 + 0.12 
2459175.57096 —24.88 + 0.13 
2459186.54032 —24.79 + 0.13 
2459257.02353 —25.33+0.11 
2459259.01527 —25.24 + 0.11 
2459269.01563 —25.43 + 0.11 
2459270.01332 -25.38 + 0.11 
2459292.94630 -25.42 + 0.11 
2459303.92858 -25.40 + 0.10 
2459313.90858 -25.44 + 0.10 
2459323.87062 -25.41 + 0.10 
2459333.96354 -25.36 + 0.10 
2459416.81508 -25.03 + 0.10 
2459434.75833 —24.84 + 0.11 
2459446.61521 —24.73 + 0.14 
2459495.57082 —24.81 + 0.13 
2459524.54114 —24.84 + 0.12 
2459536.53968 —24.88 + 0.12 
2459546.53877 —24.85 + 0.12 
2459562.54505 -24.76 +0.15 
2459629.01582 -25.27 + 0.10 


With a periodicity so close to the terrestrial year, instrumen- 
tal effects by seasonal influences must rightfully be suspected, 
and the spectrograph temperature is the first issue that comes to 
mind. Nevertheless, the variations of the spectrograph tempera- 
ture at each take of a spectrum (see Fig.|2) show no periodicity 
anywhere near a yearly cycle. Furthermore, it is important to 
note that the larger temperature variations in the first year of our 
observations was stopped by the installation of a second stage 
air conditioning of the spectrograph and its encapsulation, sig- 
nificantly improving the temperature stability for the latter years. 
There is no noticeable change or scatter of the RV measurements 
between those in the first and those in the remaining periods ob- 
served. This is additional proof that the RV recordings are real 
and not influenced by instrumental or seasonal effects. Finally, 
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Fig. 2. Temperature of the spectrograph during the TIGRE observations 
of Albireo. 


Table 2. Stellar parameters of Albireo Aa determined with an improved 


method. Values of (2021) are shown for comparison. 


Parameter value Drimmel et al. (2021) 
Teff 4358+ 10K 4382.7 +2.1 K 
logg 1.68 + 0.03 0.93 + 0.01 
[M/H] -0.08 + 0.01 0.02 
[æ/Fe] -0.10 + 0.07 0.08 
Vmic 1.85 + 0.01 km s7! 1.57 km s7! 
Viae 2.55 + 0.29 km s7! 5.22 km s~! 
vasini 4.45+0.20kms™! 8.34 = 0.4 km s”! 


we point out that in the same time interval, and with the same 
instrumentation and reduction procedures, we measured RVs for 
about 50 other stars — and none of them showed any periodicity 
close to one year. 


2.1. Stellar parameters of Albireo Aa 


We determined the stellar parameters of Albireo Aa with our 
improved method with iSpec as described in detail in 
E oR, We used the stellar parameters of 
(2021) as input parameters and performed a set of 27 fit 
calculations using the new method. This improved method more 
consistently determines the values for the surface gravity log g, 
and for the rotational and turbulence velocities. These values are 
required for our analysis. 

We present the results of our spectral analysis in Table[2] The 
effective temperature Ter has hardly changed in comparison to 
(2021). The value for the surface gravity (log g = 
1.7) is now consistent with the value calculated based on the 
parallax. The metallicity [M/H] is slightly lower. The primary 
goal was to determine the rotational velocity, which has a value 
Of va Sin i = 4.45 km s~!. Using the value for log g and assuming 
a mass of 5.2 Mo, we obtained a radius of 55 Ro for Albireo Aa. 


3. The orbital parameters of Albireo Ad 


We analyzed the RV measurements of Albireo Aa using the Ra- 


dial Velocity modeling toolkit (Fulton et al.[2018| RadVel) (ver- 


sion 1.4.9) that models RV data using the method of Bayesian in- 
ference. This toolkit fits Keplerian orbits to observed RV curves 


Table 3. Orbital parameters of Albireo Ad 


Parameter value 

P 371.2 + 5.6 days 
K 0.341 + 0.025 km s7! 
e 0.062 + 0.057 
w K 
Ta 2458867 + 18 JD 

Vrad —25.158 + 0.019 km s7! 

dydt 0.000172 + 0.000057 km s~! day! 


and determines the six orbital parameters, which are the or- 
bital period P, the time of inferior conjunction T,, eccentricity 
e, semi-amplitude K, the argument of the periapsis of the star’s 
orbit w, and the RV of the system v,aq. To obtain an estimation 
of the uncertainties, RadVel contains the Markov-Chain Monte 
Carlo (MCMC) package of Foreman-Mackey et al.|(2013). 
Because Albireo Aa is a member of a binary system with Al- 
bireo Ac, which has an orbital period of about 122 years 
[mel et al.|/2021), there is a slow change in the measured RV 


during the time of our observation campaign. The orbit deter- 
mined by (2021) has large uncertainties, espe- 
cially during the period of our observations (2018 to 2021). In 
their Figure 6, the possible orbits may differ in up to 3 km s~! 
in the RVs around the year 2020. This is larger than the semi- 
amplitude of the orbit. Therefore, we decided to fit the effect of 
an RV change induced by Albireo Ac using the dvdt parameter 
of RadVel, which fits a linear change to the observed RV curve. 

In Table |3} we list the orbital parameters determined with 
RadVel. The errors were obtained with an MCMC run using 100 
walkers and 20000 steps. The orbital period of the discovered 
companion Albireo Ad is P ~ 371 days. The semi-amplitude 
with K = 0.341 km s7! is quite small, but still above the detec- 
tion limit of 0.1 km s~! of HEROS/TIGRE (Mittag et al.[2018). 
The orbit is within the error circular (e = 0.062 + 0.057) and, 
therefore, we could not determine the argument of the periapsis 
of the star’s orbit w. It is important to note that the determined 
RV of the system (vrag = —25.158 + 0.019 km s7!) is not the sys- 
temic RV of the Albireo Aa, Ad system because it includes the 
effect of Albireo Ac. 

Considering the orbit of Albireo Ac with a period of 122 
years and a semi-amplitude of 2.91 km s7!, one obtains an aver- 
age change of 0.00026 km s~! per day, which is consistent with 
the fitted value of dvdt. The general trend of an increasing RV 
is also correct because the RV curve of Albireo Aa, Ac passed 
through the minimum a few years ago (around 2010). It is im- 
portant to keep in mind that the best fit RV curve of [Drimmel]| 
is not very well defined. 

We present the complete RV curve and the RadVel fit of Al- 
bireo Aa in Fig.|3] Subplot a shows the complete data set of the 
TIGRE RV measurements. The solid line represents the Rad Vel 
fit. The residuals are shown in subplot b. Here, the solid line 
includes the linear trend of the fitted dvdt parameter. The phase- 
folded RV curve is demonstrated in subplot c. The part around 
phase 0.0 has no RV measurements because the period of the sig- 
nal is very close to one year, and Albireo Aa was not observable 
during that phase. 

There exist published measurements of RVs of Albireo Aa 
for over a decade in the literature (Drimmel et al.|[2021). The 
oldest measurements obviously have a very large uncertainty. We 
took a closer look into the 14 RV measurements obtained with 
CORAVEL on the 1-m Swiss telescope at the Haute-Provence 


Observatory (Famaey et al.|2005). Unfortunately, the CORAVEL 
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Fig. 3. RV curve of Albireo Aa with the fit curve obtained with RadVel. 


RV data still have a scatter that is too large to detect any clear 
RV variation in the range of the amplitude of Albireo Ad. In Fig. 
we present the CORAVEL data. We corrected the data for a 
linear trend caused by Albireo Ac. The color-coding according 
to the observation epoch clearly shows that there is no other trend 
in the data. Thus, we could not use the CORAVEL data to further 
constrain the period of the Albireo Ad orbit. 


4. The mass of Albireo Ad 


We calculated the mass function f of the Albireo Aa, Ad system 

using the known formula 

f= Mi, sin? i 7 PK? a) 
(Maa + Maa) 2nG 

where G represents the gravitational constant, Ma, and M aq are 

the masses of Albireo Aa and Ad, respectively, and i is the incli- 


a- e’, 


Article number, page 4 of 6 


Year 
2021 


2022 


9000 
2450000 


9200 9400 9600 


= 371.54 days 
0.34 km s71 


Phase 


nation of the orbit. The value f = 0.00000152 + 0.00000036 Mo 
is quite small, indicating a small mass for the companion. We 
assumed a value of My, = 5.2 + 0.1 Mo for the stellar mass 
of Albireo Aa (Drimmel et al.|[2021). Solving the above equa- 
tion, we calculated the minimal mass of Albireo Ad to be 
Maasini = 0.0345 + 0.0026 Mọ. This corresponds to about 36 
Jupiter masses (M7). We assumed that the inclinations i of the 
Albireo Ac and the Albireo Ad orbits are the same. 
(2017) found a strong tendency of orbit alignment in triple 
stars when the members are within 50 AU. With the inclination 
i = 156.15° (Drimmel et al.2021), the mass of Albireo Ad is 
Maa = 0.085 + 0.007 Mo. Using the Kepler’s law equations, 
the maximal orbital separation of the two stars is about 1.9 AU. 
Assuming a distance of approximately 120 pc to Albireo A, the 
angular separation in the sky is 0.0157 arcsec. 
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Fig. 4. RV measurements of CORAVEL. Vertical axis: CORAVELe RV 
after subtraction of a linear trend. Horizontal axis: Epochs phased with 
370 days. Color: Epoch (in units of 370 days), to make any remaining 
trends visible. 


found a signal in their analysis of 
proper motion anomaly in Gaia DR2 data. They determined a 
secondary mass of m} = 416.82*!3%78M,. This signal of an 
anomaly in the proper motion can be caused by both the dis- 
tant companion Albireo Ac or the newly discovered close com- 
panion Albireo Ad. Since the proper motion anomaly is actu- 
ally a vector, the measured effect is a combination of the in- 
fluences from both companion stars. The relation between the 


Kervella et al.| (2019) mass m} and the actual mass m is given 


by m, = m / Yr, with r being the orbital radius. For Albireo Ac, 
we have r = 48 AU and m = 2.7 Mọ = 2828.44 M,. The 
contribution from Albireo Ac is, therefore, m? = 408 M,, which 
coincides with the signal found by very 
well. We performed the same exercise for Albireo Ad, with 
r = 1.9 AU and m = 0.085M, = 89 M,. The contribution 
from Albireo Ad to the signal of a proper motion anomaly found 


by |Kervella et al.|(2019) is just m? = 64.6 My. In addition, as 
described in (2019), there is a smearing factor 
y, which in the case of Albireo Ad is y = 0.1. Thus, the fi- 


nal observable signal is just 6.5 Mz. This means that the signal 
found in [Kervella et al.] (2019) almost completely originates in 
the influence from Albireo Ac, and there is no contribution from 
Albireo Ad. 


5. Discussion 


We first checked if our detected star is actually the already re- 
ported star Albireo Ab and found that this is not the case be- 
cause the angular separations are too different. We found a max- 
imal angular separation of 0.0158 arcsec for Albireo Ad from 
Aa. The angular separations reported for the detections of Al- 


bireo Ab are 0.125 arcsec (Bonneau & Foy}1980) and 0.045 arc- 
sec (Prieur et al./2002). Even the smallest measurement is al- 


ready three times larger than the maximal angular separation for 
Albireo Ad. Thus, our detected star cannot be Albireo Ab. Using 
Kepler’s third law (T? œ a°), the expected period of the Al- 
bireo Ab orbit should be at least 5.2 years. We did not detect any 
signal of that period in our RV curve, but the HEROS instrument 
is probably not sensitive enough. 

The new iSpec method for the determination of stellar pa- 


rameters of|Rosas-Portilla et al.|(2021) now gives more reliable 


1.0 


0.84 


300 400 500 600 700 


Period [d] 


Fig. 5. Lomb-Scargle periodogram of the HIPPARCOS photometry of 
Albireo A. There are no periodic variations. 
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values for the rotational velocities of stars. The determined rota- 
tional velocities have a factor of sini because of the inclination 
of the rotational axis, which is unknown. Based on the rotational 
velocity determined with iSpec of vrot = 4.45 km s7}, we calcu- 
lated the rotation period assuming a radius of 55 Ro. The result 
is a rotation period of about 618 days. This is the maximum of 
the rotation period, and depending on the inclination the rotation 
period can be shorter. 

Another way to check if we somehow detected the signal of 
stellar rotation instead of the one of the orbit of a companion star 
is by stellar activity. We checked for periodic variations in the 
common activity indicators. The TIGRE telescope has a data re- 
duction pipeline that determines the S-index (Mittag et al.[2016), 
which is a measure for the magnetic activity of the Can H&K 
lines. We found a small peak (~ 0.4) in the Lomb-Scargle pe- 
riodogram at 370 days. However, this part of the spectrum is 
completely dominated by the flux of Albireo Ac. The routine 
that measures the S-index obviously corrects for the RV of Al- 
bireo Aa. This introduces variations into the S-index, because 
the flux of Ac is falsely shifted in wavelength and, therefore, in 
the window that is used to determine the S-index. Because the 
RV correction has the variation of the orbital period, it obviously 
causes this false signal. 

Variations of the RV in stars caused by stellar activity (star 
spots) are a widely known problem for the detection of exoplan- 
ets. Several studies have been performed to quantify this effect 
[Desort et al.|2007}. As in the case of 
a companion, the signal is a sine curve, but there are usually 
jumps (phase shifts) in the signal because of the appearance and 
disappearance of the star spots or active regions. We did not de- 
tect any phase shifts during over three years of RV observations. 
The expected variations caused by stellar activity have normall 
peak-to-peak amplitudes of at most 0.1 km s7! 
[1997). This is below our detected RV signal with a peak-to-peak 
amplitude of 0.68 km s~!. Thus, the detected signal of RV varia- 
tions is caused by the orbital movement of a companion star and 
cannot come from stellar activity on Albireo Aa. 

We may comfortably exclude the possibility of a pulsation 
or oscillation for the RV signal. The former possibility can be 
excluded by the absence of a matching photometric signal in 
the HIPPARCOS data (van Leeuwen|2007). If there were any 
pulsations of Albireo Aa, then a brightness change should be 
observed. As shown in Fig.|5} the periodogram of the HIPPAR- 
COS photometry data has no peak at all. This can also be seen in 
Fig. [6| where the phase-folded photometry of Albireo A is pre- 
sented. The idea of an oscillation is excluded by the expected 
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Fig. 6. HIPPARCOS photometry data of Albireo A folded with the 
371 days period show no variations. 


period being much shorter — after all, Albireo Aa is a much 


more compact giant than a Mira star. According to|Kjeldsen & 
Bedding}(1995), the frequency of the maximum-amplitude oscil- 


lation, Vmax, Of a star (which is the well-known 5-minute mode 
of the Sun) mainly depends on gravity g and effective temper- 
ature Ter by being proportional to g/ VTer. Using the physical 
parameters of Albireo Aa listed in Table [2] Vmax Of Albireo Aa 
is exactly 500 times smaller than the 5-minute oscillation of the 
Sun, meaning a period of 1.74 days. 


The orbital period of 371 days is very close to one year. Be- 
cause Albireo A is at least a triple system, there must be a large 
difference between the orbits of Ac and Ad, so that the system 
is gravitationally stable. The existence of the star Albireo Ab 
would actually disturb this configuration, making the system un- 
stable. 


(2021) found evidence that Albireo Ac itself 


is a binary system. They determined the mass ratio of Albireo Aa 
and Ac, and it indicates a "missing mass" for Albireo Ac. In this 
work, we discovered a new star in the Albireo A system that 
obviously changes the mass distribution. However, the additional 
mass is very small (0.085 Mọ) and has no effect on the mass 
problem of the Albireo Aa, Ac system. 


In passing, we note that — depending on the unknown ori- 
entation of the Albireo Ad orbit, that is to say depending on its 
projection onto the celestial sphere — the 371-day period does 
modify the astrometrically measured parallax of Albireo A. If 
the inclination of the orbit is indeed about 156 degrees and if the 
semi-major axis is thus about 16 mas, then the radius of the Aa 
orbit around the common center of mass is about 0.25 mas. If the 
inclination is farther from 90 or 270 degrees, that is if the actual 
mass of Ad is higher than 0.085 Mo, the astrometric orbit of Aa 
is correspondingly larger. The parallaxes of Albireo A and B, as 
published in Gaia eDR3, are 8.98+0.45 mas and 8.19+0.08 mas, 
respectively. The difference of 0.79 mas is a bit less than twice 
the combined uncertainty. The astrometric effect of the newly 
discovered component Ad might possibly explain a significant 
part of this difference. 
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6. Conclusions 


The famous bright multiple stellar system Albireo remains an 
object of interest, to which it is worth dedicating further observa- 
tions and studies. In our analysis of the RV curve of Albireo Aa, 
we detected a clear signal of a new member of the system. Al- 
bireo Ad is a star with a mass of 0.085 Mo orbiting the giant star 
Albireo Aa with a period of P = 371 days. The orbit is close to 
circular and has a small semi-amplitude of 0.34 kms. 

We conclude that the Albireo system has a clear hierarchical 
structure. There is the possible gravitational connection between 
Albireo A and B. Albireo A is a multiple stellar system. The 


companion Albireo Ac is probably a binary star (Drimmel et al. 
2021). In addition, the new member Albireo Ad has a close orbit 


to Albireo Aa. With the discovery of the new star Albireo Ad, 
the existence of Albireo Ab is very unlikely. 


Acknowledgements. This research has been made possible by the CONACyT- 
DFG bilateral grant No. 278156. We thank the University of Guanajuato for 
the grants for the projects 036/2021 and 105/2021 of the Convocatoria Insti- 
tucional de Investigacion Cientifica 2021 and 2022. This work has made use of 
data from the European Space Agency (ESA) mission Gaia 


cosmos.esa.int/gaia), processed by the Gaia Data Processing and Anal- 
ysis Consortium (DPAC, https: //www.cosmos.esa.int/web/gaia/dpac/ 
. Funding for the DPAC has been provided by national institutions, 
in particular the institutions participating in the Gaia Multilateral Agreement. 
This research has made use of the VizieR catalogue access tool, CDS, Stras- 


bourg, France (Ochsenbein et al.|2000). 


References 


Bastian, U. & Anton, R. 2018, A&A, 620, L2 

Bonneau, D. & Foy, R. 1980, A&A, 86, 295 

Clerke, A. M. 1899, The Observatory, 22, 387 

Desort, M., Lagrange, A. M., Galland, F., Udry, S., & Mayor, M. 2007, A&A, 
473, 983 

Drimmel, R., Bucciarelli, B., & Inno, L. 2019, Research Notes of the American 
Astronomical Society, 3, 79 

Drimmel, R., Sozzetti, A., Schröder, K.-P., et al. 2021, MNRAS, 502, 328 

Famaey, B., Jorissen, A., Luri, X., et al. 2005, A&A, 430, 165 

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 
306 

Fulton, B. J., Petigura, E. A., Blunt, S., & Sinukoff, E. 2018, PASP, 130, 044504 

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2021, A&A, 650, C3 

Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, Al 

Griffin, R. F. 1999, JRASC, 93, 208 

Hartkopf, W. I. 1999, IAU Commission on Double Stars, 139, 2 

Hass, J., Phung, K., Carro, J., et al. 2016, Journal of Double Star Observations, 
12, 204 

Jack, D., Hernandez Huerta, M. A., Rosas Portilla, F. D., & Schréder, K.-P. 2022, 
Astronomische Nachrichten, accepted 

Jack, D., Hernández Huerta, M. A., & Schröder, K.-P. 2020, Astronomische 
Nachrichten, 341, 616 

Jack, D., Schröder, K.-P., & Bastian, U. 2018, Research Notes of the American 
Astronomical Society, 2, 225 

Kervella, P., Arenou, F., Mignard, F., & Thévenin, F. 2019, A&A, 623, A72 

Kjeldsen, H. & Bedding, T. R. 1995, A&A, 293, 87 

Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021, A&A, 649, A2 

Markowitz, A. H. 1969, PhD thesis, THE OHIO STATE UNIVERSITY. 

Mason, B. D., Hartkopf, W. I., & Hurowitz, H. M. 2013, AJ, 146, 56 

Maury, A. C. & Pickering, E. C. 1897, Annals of Harvard College Observatory, 
28, 1 

Mittag, M., Hempelmann, A., Fuhrmeister, B., Czesla, S., & Schmitt, J. H. M. M. 
2018, AN, 339, 53 

Mittag, M., Schröder, K. P., Hempelmann, A., González-Pérez, J. N., & Schmitt, 
J. H. M. M. 2016, A&A, 591, A89 

Ochsenbein, F., Bauer, P., & Marcout, J. 2000, Astronomy and Astrophysics Sup- 
plement Series, 143, 23 

Parsons, S. B. & Ake, T. B. 1998, ApJS, 119, 83 

Prieur, J. L., Koechlin, L., Ginestet, N., et al. 2002, ApJS, 142, 95 

Roberts, L. C. & Mason, B. D. 2018, MNRAS, 473, 4497 

Rosas-Portilla, F. D., Schröder, K.-P., & Jack, D. 2021, MNRAS, submitted 

Saar, S. H. & Donahue, R. A. 1997, ApJ, 485, 319 

Scardia, M. 2019, IAU Commission G1, 198, 3 

Scardia, M., Prieur, J. L., Pansecchi, L., et al. 2007, MNRAS, 374, 965 

Schmitt, J. H. M. M., Schröder, K.-P., Rauw, G., et al. 2014, AN, 335, 787 

Tokovinin, A. 2017, ApJ, 844, 103 

van Leeuwen, F. 2007, A&A, 474, 653 


